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Abstract 

Correlating expressions for two-phase flow breakup 
of liquid nitrogen, LN 2 , jets in sonic velocity nitrogen 
gasflows were obtained for an atomizing- gas temperature 
range of 111 to 442 K. The correlations were based on 
characteristic dropsize measurements obtained with a 
scattered-light scanner. The effect of droplet vaporiza- 
tion on the measurements of the volume-median drop- 
size, D y 5 , was calculated by using previously determined 
heat and momentum transfer expressions for droplets 
evaporating in high-velocity gasflow. Finally, the drop- 
size of the originally unvaporized spray, D v 5c , was cal- 
culated, normalized with respect to jet diameter, D Q , 
and correlated with atomizing-gas flowrate and tempera- 
ture, according to the following expression: 
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where WeRe = P^D^Vc/ • Here is liquid viscosity, 
p g and p 1 are gas and liquid densities, respectively, a is 
surface tension, V c is acoustic gas velocity and T g is 
atomizing-gas temperature, normalized with respect to 
airstream temperature, T Q = 293 K. This expression 
agrees well with atomization theory which predicts 
D y 5 ~ Vj* 33 , for liquid-jet breakup in high-velocity 
gasflow. 

Nomenclature 

A 0 fuel nozzle orifice area, cm 2 
a acceleration, cm/sec 2 

C d drag coefficient 

D 0 liquid-jet diameter, cm 
D v 5 volume median drop diameter, cm 
k correlation coefficient for Eq. (l) 

k' correlation coefficient for Eq. (6) 

k" correlation coefficient for Eq. (7) 

Nu heat-transfer Nusselt number, based on D y 5e 
n exponent for Eq. (1) 


Re Reynolds number based on D y 5e 
t vaporization time, sec 

V c acoustic velocity, cm/sec 

W weight flow of fluid, g/sec 
We Weber number based on D y 5e 

x axial downstream spray sampling distance 

a surface tension relative to air, dynes/cm 
\l absolute viscosity, g/cm sec 

p fluid density, g/cm 3 

Subscripts 
c calculated 

d droplet 

e experimental 

f free-stream 

g gaseous nitrogen, GN 2 

1 liquid nitrogen, LN 2 

o orifice 

Introduction 

An experimental investigation of cryogenic liquid- 
jet breakup in high-velocity gasflow wets conducted to 
determine the effect of atomizing-gas temperature, T g , 
on characteristic dropsize, D y 5 , of liquid-nitrogen 
sprays. Very little data are available in the spray 
literature that show gas-temperature effects on atomi- 
zation. However, the effect of gas velocity, V g , on 
liquid-jet breakup in gas streams has been studied by 
numerous investigators. 1-7 Their results are summarized 
in Table I. Water-jet breakup results described in Ref. 1 
show that good agreement can be obtained with atomi- 
zation theory, when dropsize measurements are made 
close to the atomizer orifice. However, a marked effect 
of droplet vaporization on dropsize measurements of 
water sprays did occur, as reported in Ref. 1, when the 
sampling distance downstream of the atomizer was in- 
creased from 2.2 to 6.7 cm. In the case of liquid-nitrogen 
sprays, good agreement with atomization theory may 
not occur, since highly volatile sprays can quickly 
become partially vaporized. 
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In a previous investigation of disintegrating LN 2 
jets, 2 dropsize measurements were made in the presence 
of relatively high thermal gradients. The atomizing gas, 
GN 2 , was at room temperature, 293 K, whereas LN 2 
droplet temperatures were near the boiling point of LN 2 , 
77 K. Thus, heat transfer across the gas film had a 
driving potential, AT, of 216 K. This is considerably 
higher than values of AT in the order of 10 to 15 K, 
which were encountered in previous water spray studies 
reported in Ref. 1. 

Although the original characteristic dropsize 
initially formed at the atomizer orifice was not measured 
directly, in the present study, it was possible to deter- 
mine the initial value of D y 5c by calculating the 
amount it had changed due to evaporation. Droplet ac- 
celeration and vaporization rates were calculated from 
heat transfer and drag coefficient data given in Ref. 8. 
These data had previously been obtained for drops accel- 
erating and vaporizing in high-velocity gasflows. Such 
data are difficult to obtain for drops that are micro- 
scopic in size and attain high velocity in a short distance 
of travel. By using heat transfer and drag coefficients 
given in Ref. 8, values of AD y 5 produced by droplet 
vaporization could be calculated and used to determine 
original values of D y 5c , that existed before vaporization 
and occurred. 

To determine the effect of atomizing- gas tempera- 
ture and gas mass-flux on spray dropsize produced by 
liquid-jet breakup in high-velocity gasflow, the charac- 
teristic dropsize D y 5e was measured with a scattered- 
light scanner developed at NASA Lewis Research Center 
by Buchele. 9 Sprays were sampled with the laser beam 
center line positioned at a distance of 1.2 cm down- 
stream of the fuel nozzle orifice to minimize the loss of 
small droplets due to vaporization. Volume mean diam- 
eter, D y 5e varied from 3 to 30 fim and measurements of 
D v.5e were made at atomizing- gas temperatures of 111, 
293, and 422 K, respectively. 

Apparatus and Procedure 

A two-fluid nozzle was used with assist nitrogen 
gasflow, GN 2 , to breakup liquid nitrogen, LN 2 -jets, as 
shown in Fig. 1. It was mounted at the center line of the 
24-cm diameter duct and operated over pressure ranges 
of 0.2 to 1.0 MPa for both LN 2 and GN 2 . LN 2 sprays 
were injected downstream into the airflow, just up- 
stream of the duct exit, and sampled at a distance of 
1.2 cm downstream from the atomizer orifice to the cen- 
ter line of the 4.4 X 1.9 cm laser beam. The two-fluid 
nozzle was fabricated according to the diagram illustra- 
ted in Fig. 2. LN 2 at a temperature of 77 K was axially 
injected into the airstream by gradually opening the 
control valve until the desired flowrate of 51 g/sec was 


obtained as indicated by a turbine flowmeter. The atom- 
izing gas was then turned on and weight flowrate was 
measured with a 0.51-cm diameter sharp-edge orifice. 
After the air, GN 2 and LN 2 flowrates were set, the 
volume median diameter, D y 5e , was measured with the 
scattered-light scanner. 

The optical system of the scattered-light scanner 
shown in Fig. 3 consisted of a laser beam expander with 
spatial filter, rotating scanning-slit and a detector. The 
instrument measures scattered light as a function of 
scattering angle by repeatedly sweeping a variable- 
length slit in the focal plane of the collecting lens. The 
data obtained is scattered-light energy as a function of 
the scattering angle relative to the laser beam axis. This 
method of particle size measurement is similar to that 
described in Ref. 10. Measurements of scattered-light 
energy normalized by the maximum energy were plotted 
against scattering angle and used to determine volume- 
median diameter, D y 5e , as described in Ref. 11. Also, it 
should be noted that the size-distribution dispersion can 
also be determined from this plot. Also, this method of 
determining characteristic dropsize and dispersion of 
dropsize can be used independent of particle size distri- 
bution function, according to Buchele. 11 For a typical 
measurement, the scan is repeated 60 times per second 
to average out any temporal variations in the energy 
curve. 

Spray pattern effects were minimized by measuring 
D y 5e for the entire cloud of droplets. The instrument 
was calibrated with five sets of monosized polystyrene 
spheres having diameters of 8, 12, 25, 50, and 100 pm. 
Since the sprays were sampled very close to the atomizer 
orifice, they contained a relatively high number-density 
of very small droplets. As a result, the light-scattering 
measurements required correction for multiple scattering 
as described in Ref. 12. Also, dropsize measurements 
were corrected to include Mie scattering theory when 
very small drop diameters, i.e., 10 /*m, were measured. 
Reproducibility tests showed that experimental measure- 
ments agreed within ±5 percent. 

It was necessary to correct experimental meas- 
urements by taking background readings when the 
atomizing-gas, GN 2 temperature was well above or 
below airstream temperature. This was due to high gas- 
density gradients being present when atomizing-gas tem- 
peratures were relatively high, 422 K, or low, 111 K, as 
compared with the surrounding airstream temperature 
of 293 K. Only a small correction was needed at 
T g = 293 K, since gas density gradients were close to 
zero. Temperature gradients for nitrogen vapor films 
surrounding the droplets were assumed to have negligi- 
ble effect on measurements, since this was found to be 
the case when water sprays were studied in Ref. 15. 


2 


Experimental Results 

To obtain dropsize measurements, the entire LN 2 - 
spray cross section was sampled and the laser beam 
center line was located at a distance of 1.2 cm down- 
stream of the fuel-nozzle orifice, as shown in Fig. 1. 
Droplets traveled a distance of 2.0 cm in passing 
through the scattered-light scanner laser beam. However, 
some of the very small and highly volatile LN 2 droplets 
were completely vaporized before they could travel 
through the laser beam. As a result, experimental values 
of D y 5e were obtained for partially vaporized sprays. 
Thus, it was necessary to calculate the change in 
dropsize, AD y 5 , in order to estimate the initially 
unvaporized spray dropsize, D y 5c , that was formed at 
the fuel-nozzle orifice. Values of D y 5c , were then 
correlated with atomizing-gas flowrate, W g . Such 
correlations are needed for modeling spray vaporization 
and combustion processes. 

Effect of GN ? Flowrate on D y 5e 

Measurements of D y 5e were made with the 
scattered-light scanner and plotted against GN 2 flow- 
rate, W g , as shown in Fig. 4. Since high-velocity 
atomizing-gas flowrates were used, cryogenic liquid-jet 
breakup occurred primarily in the regime of aerody- 
namic stripping. No indication of secondary breakup of 
droplets was observed since low atomizing-gas velocities 
were not used. Thus, the low gas-velocity regime of cap- 
illary wave breakup of liquid jets was not investigated. 

From the plot shown in Fig. 4, reciprocal D y 5e was 
correlated with atomizing-gas flowrate, W g , and the 
following expression was obtained: 


did appear to agree with theory, the proportionality 
constant k was too low to adequately characterize the 
initial unvaporized spray. Thus, the study in Ref. 9 did 
not take into account the effect of small droplets vapor- 
izing completely before they could pass through the laser 
beam. 

Acceleration of LN ? Droplets 

The effect of droplet vaporization rate on experi- 
mental values of D y 5e was determined by calculating 
vaporization time, t, as based on droplet velocity V d , for 
D y 5e . Time, t, was calculated over a distance of 2.2 cm, 
i.e., the distance from nozzle orifice to the downstream 
edge of the laser beam, as shown in Fig. S. 

In order to determine volume-median drop veloc- 
ity, V d , the acceleration, a, of LN 2 droplets was calcu- 
lated from the following momentum balance as given in 
Ref. 8: 

m d a = lp g A d (V g - V d ) 2 C d (2) 

2 

where m d and A d are mass and area of dropsize D y 5e , 
respectively, i.e., m d = 5e /6, ^ ^e drag 

coefficient based on characteristic length, D y 5e . Rewrit- 
ing this expression, in terms of change in drop-velocity 
squared, AV d , over the distance of travel, Ax, the 
following relationship is obtained: 

iv 2 , V,(V, - V d )’ r (3| 

A* Vi Dv.b. “ 


Values of the proportionality constant k and exponent 
n are given in Table II. At an atomizing-gas tempera- 
ture of 293 K, the following expression was obtained: 
D v * 5e = SOI where D y 1 5e and W g are expressed 

as cm -1 and g/sec, respectively. 

The exponent 1.11 for W is considerably less than 
1.33 as predicted by theory for liquid-jet breakup in 
high-velocity gasflow. This discrepancy can be attrib- 
uted to a loss of small vaporizing LN 2 drops before 
spray measurements could be made with the scattered- 
light scanner. In the present study, results agree better 
with atomization theory than those reported in Ref. 9. 
This is due to the allowance made in the present study 
for the effect of droplet vaporization on dropsize meas- 
urements of highly volatile sprays. This effect was not 
accounted for in Ref. 9 and although the dropsize data 


where C d = 27 Re 0,84 , as given in Ref. 8, and Re is 
based on the characteristic dropsize, D v 5c . 

The deceleration of atomizing gaseous nitrogen jets 
into a surrounding low velocity airflow was determined 
as follows. At the nozzle orifice, gas velocity, V , was 
equal to the acoustic velocity, V c , of gaseous nitrogen. 
Values of V g used to solve Eq. (3), were calculated at 
downstream distances of x = 0.5, and 10 cm, respec- 
tively and plotted in Fig. 5. Calculated values of V g are 
based on data given in Ref. 13 and plotted in Fig. 5 for 
comparison. The percent deceleration of the atomizing 
nitrogen gas is assumed to be approximately the same in 
both cases, since the two-fluid nozzles used in Ref. 13 
and the present study were very similar in design. 

To determine the acceleration of LN 2 droplets 
characterized by D y 5e , values of V d were calculated by 
numerically integrating Eq. (3) and plotting V d against 
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downstream distance, x, as shown in Fig. 6, for three 
atomizing- g £is temperatures. LN 2 droplet vaporization 
time, At, was calculated from this plot by means of the 
expression At = Ax/V d . Calculated values of At for a 
given distance Ax are given in Table ID along with cal- 
culated Reynolds numbers averaged over the incremental 
distance Ax and values of D y 5c . Gas and liquid trans- 
port properties used in calculating vaporization times 
are given in Table IV. 


Cryogenic Spray Vaporization Rates 

Vaporization rates of LN 2 sprays characterized by 
D v 5 e were calculated by using the following heat- 
balance equation: dm d /dt = hA AT/H t , where h is the 
heat-transfer coefficient, A is droplet surface-area based 
on D v 6e , AT = T g - T d , and H* = 1^ + C p AT. Here 
Hy is the latent heat of vaporization of LN 2 and C p is 
the specific heat of GN 2 . This expression may be rewrit- 
ten as follows to obtain vaporization rate in terms of 
changes in droplet surface-area with time: 

AD v.5 = a TNu (4) 

At p 1 H t 

where k g and p 1 are gas thermal conductivity and liquid 
density, respectively. In previous fuel droplet studies 
reported in Ref. 8, a high-speed droplet tracking camera 
was used to determine vaporization rates of fuels such 
as n-octane, jet-A, and numerous other liquids includ- 
ing water, benzene, acetone, and carbon tetrachloride. 
It was found that: Nu = 2 + 0.30S Re 0 ' 6 , where Re 
= p g D y 5e A V//i g) and AV is the average velocity dif- 
ference over an incremental distance Ax. Vaporization 
rate calculations were based on the characteristic drop 
diameter D y 5c . GN 2 viscosity and thermal conductivity 
were evaluated at the average gas-film temperature, i.e., 
Tf = 1/2 (T + TJ. LN 2 droplet surface temperatures 
were assumed to be near the boiling point, 77 K, when 
droplet sprays were being accelerated and partially 
vaporized. The latent heat of vaporization of LN 2 was 
evaluated at 77 K and the specific heat of nitrogen va- 
por was evaluated at the average gas-film temperature, 
Tf- 

Experimental values of D~* 5e , calculated values of 
D” 1 5c and changes in characteristic drop diameter 
squared, AD y 5 , that occurred due to partial vaporiza- 
tion of the cryogenic sprays are given in Table V. Cal- 
culations were based on GN 2 and LN 2 flowrates of 
4.54 and 51 g/sec, respectively. Values of D y 5c were 
calculated from the expression: — AD y 5 = D y 6c — D y 5e 
and they are plotted against W g as shown in Fig. 7. 


From this plot, the following correlating expression is 


obtained: D y 5c = k c W g , which is similar to Eq. (l). 
Values of k c and n are given in Table V. For the case T 
= 293 K: 


D 1 = 285 W 1-33 
v.5c g 


( 5 ) 


Comparing this expression with Eq. (1), i.e., D“ 1 5e 
= 301 W g * 1 , shows that the proportionality constants k e 
and k c are nearly equal. Also, Eq. (5) shows that droplet 
vaporization did have considerable effect on the expo- 
nent n. The value of n given in Eq. (5) agrees well with 
atomization theory which predicts n = 1.33. Thus, the 
agreement of Eq. (5) with atomization theory indicates 
that an expresson for the unvaporized spray near the 
fuel-nozzle orifice can be calculated using heat-transfer 
and drag coefficients given in Ref. 8. 


Correlation of D y with Dimensionless Force Ratios 

The calculated volume- median diameter D y 5c was 
normalized with respect to LN 2 jet diameter, D 0 , and as 
shown in Fig. 8, is plotted against the product of We, 
Re, and Pg/Pp i.e., the Weber number, Reynolds 
number, and gas-to-liquid density ratio, respectively. 
From this plot, the following dimensionless expression 
was derived: 



( 6 ) 


where WeRe is the ratio of aerodynamic to liquid-jet 
surface forces, i.e., liquid viscosity and surface tension. 
From the three plots shown in Fig. 8, it is evident that 
D y 5c is a function of atomizing-gas temperature, i.e., 
D Q /D y 5c = f(T g ). Thus, k' is assumed to be a function 
of atomizing-gas temperature normalized with respect to 
Tf, i.e., k' ~ T /Tf, where Tf is the free-airstream tem- 
perature, 293 K, of the low velocity airflow in the test 
section as shown in Fig. 1. 

From a log-log plot of k' against T g /T f , it was 
found that k' ~ (T g /Tf ) 1,25 and as a result, the follow- 
ing correlating expression can be written: 
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v.5c 
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0.44 (r. 


1.25 


( 7 ) 


This expression is plotted in Fig. 9 which shows that k£ 
= 9.0. Thus, Eq. (7) may be rewritten as follows: 
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where p g V c = W g /A Q . Here it is evident that the ther- 
modynamic effect of normalized atomizing-gas tempera- 


ture, T g /Tf on the reciprocal volume-median diameter 
is nearly as great as that of the atomizing-gas mass flux 
which is also a function of gas temperature, i.e., 
As a result, Eq. (8) shows that D” 1 


pV’~ T -0 - 5 . 

b 58 g 

~ rpu.oo ^ similar effect of atomizing-gas temperature 
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on D y 5 was obtained for water sprays, as reported in 
Ref. 15. 


From experimental dropsize measurements of 
partially vaporized liquid-nitrogen sprays, it was found 
that increasing the atomizing-gas temperature gave a 
marked increase in the surface area per unit volume of 
liquid-nitrogen sprays. A result that would be very 
beneficial in producing very rapid and efficient fuel- 
spray combustion in gas turbine and rocket combustors. 


Concluding Remarks 

Computations in the present study were based on 
numerical integration of momentum and heat-transfer 
expressions that had been developed in previous droplet 
studies, at NASA Lewis. As a result, the effect of droplet 
vaporization on dropsize measurements was determined. 
Without this source of knowledge, on calculating rates 
of heat and momentum transfer to vaporizing cryogenic 
drops, it would be almost impossible to determine the 
characteristic dropsize of an initially unvaporized 
cryogenic spray produced at the orifice of a two-fluid 
fuel nozzle. Thus with the computational method used 
in this study, it was possible to determine the effect of 
atomizing-gas temperature on cryogenic spray dropsize, 
once the effect of droplet vaporization on dropsize 
measurements had been determined. Also, it was found 
that the effect of atomizing-gas mass-flux, p g V c , on 
volume-median drop diameter agreed well with atomiza- 
tion theory for liquid-jet breakup in high-velocity 
gasflows. The final correlating expression derived in this 
study (Eq. (8)), can be readily applied to analytical 
modeling of fuel spray vaporization and combustion in 
gas-turbine and rocket combustors, within the range of 
variables investigated in this study. 
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TABLE I. — ATOMIZING-GAS VELOCITY EXPO- 
NENT, n, FOR HIGH-VELOCITY GASFLOW 


BREAKUP OF LIQUID JETS 


Source 

Exponent, 

n 

Adelberg, Theory 14 

1.33 

Present study, x = 2.2 cm 

1.33 

Kim and Marshall 3 

1.14 

Lorenzetto and Lefebvre 4 

1.00 

Nukiyama and Tanasawa, 5 x = 5 to 25 cm 

1.00 

Weiss and Worsham 6 

1.33 

Wolf and Anderson 7 

1.33 


TABLE H.— COEFFICIENT k 
AND EXPONENT n FOR 


EQ. (1)» 


Atomizing- gas 
temperature, 

T «’ 

K 

k 

n 

111 

82 

1.22 

293 

301 

1.11 

422 

367 

1.08 


*D 


— l 
v.5e 


= k t W 


n 

g 


TABLE III.— VAPORIZATION TIME, t, 

FOR D~g 4 AT W g = 4.54 g/sec AND 
Ax = 2.2 cm 


Atomizing- gas 

D -1 

AtxlO 4 , 

temperature, 

v.5e’ 

sec 

T «- 

cm"” 1 


K 



111 

540 

4.40 

293 

1600 

1.44 

422 

1900 

.99 
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TABLE IV— GN 2 AND LN 2 TRANSPORT PROPERTIES AT 

W, = 4.54 g/gec AND T„ = 293 K 

Nitrogen gag 


V c , cm /gee 3.43X10 4 

P t , g/cc 3.84xl0“ 4 

g/cm gee 1.25x10 4 

'k g , eal/gee gq cm (*C/cm) 4.20x10“ 5 

LN 3 droplets 

V d , cm/sec 204 

P v g/cc 0.80 

H v , cal/g 47.8 

Cn V , cal/g *C 0.25 


'Evaluated at T f = l/2(T g - T x ). 


TABLE V.— REDUCTION IN DROPSIZE, -A D* 8 , AND 


UNVAPORIZED DROPSIZE, D v g , AT Vf = 4.54 g/gec 


Atomizing- gas 
temperature, 

T «. 

K 

-AD*. # xl0 9 , 

cm 2 

cm -1 

d ;.sc 

cm"” 1 

111 

29.7 

540 

565 

293 

16.3 

1600 

2100 

422 

13.2 

2000 

. 2950 



Figure 1 .—Apparatus and auxilliary equipment 
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Atomizing gas * 




Inside diameter 
0.584 cm 

\ 


m 










1 Nozzle orifice 

Uquid nitrogen diameter 

0.663 cm — ^ 

Figure 2.— Diagram of pneumatic two-fluid atomizer. 


Atomizer, D 0 = 0.564 cm 


Light 
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angle .13 rad 



(-« — 10 cm 
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multiplier 

Motor 

'‘—Timer 
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Figure 3. — Atmospheric pressure test section and optical path of scattered-light scanner. 



Figure 4. — Effect of atomizing-gas flowrate on reciprocal volume- 
median diameter, D" 1 ^, for partially vaporized LN 2 sprays. 


8 


Atomlzing-gas velocity, V§, (cm/sec) 



x, cm 

Figure 5. — Deceleration of atomizing -gas, GN 2 , 
downstream of fuel -nozzle orifice. 
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Atomizing -gas 
flowrate, 


D ;ie- cm_1 



10x10® Atomizing-gas 


D Y.5e* cm ' 





(c) Tg = 422 K. 


Figure 6.— Acceleration of volume-median dropsize, D v . 5 e . 
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Do/D v .5c <* <o Reciprocal volume-median diameter, D~] gc 



ure 7. — Calculated effect of Wg on initial unvaporized value of 
5 C , at fuel nozzle orifice, x = 0. 


Tg. K k 


O 422 13.7 



Figure 8.— Correlation of D^/Dy Sc with dimensionless groups. 



Figure 9. — Correlation of Dc/Dyse with dimensionless groups. 
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